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MAXIMUM  PERIODIC  WAVE  RUNUP  ON  SMOOTH  SLOPES 


Introduction 

Wave  runup  is  defined  as  the  height  above  Stillwater  level  to  which  a 
wave  will  rise  on  a  structure  or  beach.  Proper  design  of  coastal  structures 
depends  on  the  ability  to  make  reliable  estimates  of  the  maximum  runup  that 
might  occur  on  the  structure  during  a  significant  storm  or  extreme  event.  A 
method  is  presented  in  this  paper  which  provides  a  conservative  approach  for 
structures  fronted  by  flat  slopes. 

Considerable  effort  has  addressed  the  calculation  of  wave  runup  on 
structures  due  to  short  (wind  and  swell)  wave  action.  Summaries  of  studies  on 
wave  runup  for  periodic  waves  using  various  methods  of  data  presentation  were 
reported  in  Saville  (1956),  Savage  (1958),  Koh  and  LeMehaute  (1966),  Van  Dorn 
(1966),  Raichlen  and  Hammock  (1974),  Battjes  (1974),  Technical  Advisory  Com¬ 
mittee  on  Protection  Against  Inundation  (1974),  and  Stoa  (1978).  Numerous 
other  references  have  also  been  written  on  the  subject  dealing  with  irregular 
wave  runup  and/or  long  wave  runup.  The  present  paper  presents  another 
approach  to  data  presentation  for  the  case  of  short  monochromatic  waves  on 
smooth  surface  sloping  structures  with  an  intent  to  unify  various  existing 
theories  and  approaches  for  wave  runup  and  to  provide  a  reasonable  means  to 
calculate  an  upper  limit  of  runup  for  conservative  design  practice.  A  partic¬ 
ular  shortcoming  in  much  of  the  early  work  on  runup  in  the  US  as  summarized  in 
the  present  runup  design  curves  presented  in  the  US  Army  Corps  of  Engineers 
Shore  Protection  Manual  (1984)  is  that  authors'  interpretations  of  data  curves 
have  been  shown  and  the  data  itself  deleted  from  such  curves.  The  lack  of 
data  points  does  not  allow  the  engineer  to  estimate  accuracy  from  what  scatter 
might  exist.  The  reanalysis  of  periodic  wave  runup  data  provided  in  this 
paper  provides  the  raw  data  points  in  a  new  format.  An  additional  variant  in 
the  present  reanalysis  of  runup  data  is  to  provide  wave  runup  in  terms  of  wave 
height  at  the  structure  toe  depth  as  opposed  to  utilizing  the  deep  water  wave 
height.  The  advantage  of  utilizing  transformed  wave  height  rather  than  deep 
water  wave  height  is  that  wave  height  transformation  uncertainty  from  deep 
water  to  the  structure  site  becomes  a  separate  problem,  uncoupled  from  the 
problem  of  what  runup  actually  occurs  on  the  structure  due  to  a  given  wave  at 
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the  toe  of  the  structure.  Although  it  Is  recognized  that  the  type  (shape)  of 
wave  existing  at  the  site  is  important  to  the  ultimate  problem  (i.e.  the 
transformation  prior  to  the  structure  and  the  ultimate  runup  are  not  entirely 
uncoupled) ,  the  presentation  of  data  in  terms  of  wave  conditions  at  the  base 
of  the  slope  should  be  of  benefit  to  the  ultimate  user  of  runup  curves  who 
might  want  to  design  using  depth  limited  breaking  wave  height  criteria.  Since 
wave  period  is  considered  invariant  throughout  the  transformation  process, 
deepwater  wave  length  is  still  used  in  the  analysis. 

Data  Sources 


The  data  sources  for  this  runup  reanalysis  are  from  earlier  tests  at  the 
Coastal  Engineering  Research  Center  on  smooth  slope  runup.  These  data  are 
discussed  in  length  by  Saville  (1956)  and  Savage  (1958).  This  same  set  of 
data  was  used  by  Stoa  (1978)  in  an  earlier  reanalysis  of  runup.  For  purposes 
of  clarification,  a  short  discussion  of  this  data  set  follows.  Further  infor¬ 
mation  on  these  tests  can  be  found  in  Saville  (1956)  and  Savage  (1958). 

The  test  procedure  involved  placing  a  smooth  surface  plywood  test  slope 
in  the  end  of  the  wave  tank  and  propagating  periodic  waves  of  known  character¬ 
istics  toward  the  slope.  The  waves  in  each  test  were  measured  after  the  ini¬ 
tial  unsteady  wave  transients  died  down  but  prior  to  rereflection  of  waves 
from  the  wave  generator.  An  average  of  six  to  fifteen  waves  were  visually 
measured  by  reading  the  runup  on  a  scale  marked  on  the  face  of  the  slope  to 
the  nearest  hundredth  of  a  foot  in  vertical  elevation.  In  all  the  data  pre¬ 
sented  in  this  paper,  the  water  depth  was  constant  (-  0.38  meters).  Saville 
(1956)  noted  that  varying  the  water  depth  at  the  toe  of  the  structure  had  a 
negligible  effect  on  the  relative  wave  runup  when  the  water  depth  at  the  toe 
of  the  structures  was  equal  to  or  greater  than  three  times  the  deep  water  wave 
height. 

Wave  characteristics  were  determined  by  calibrating  the  wave  generator 
for  the  0.38  meter  water  depth.  The  generator  was  calibrated  by  placing  a 
wave  absorber  in  the  beach  end  of  the  tank  and  generating  a  wave  train  with 
known  and  reproducible  settings  on  the  generator.  The  average  height  of  the 
wave  train  so  generated  was  measured  with  a  parallel  wire  gage  at  2  meter 
intervals  along  the  tank  beginning  near  the  wave  generator.  Wave  heights  were 
plotted  versus  distance  along  the  tank  and  the  wave  height  value  obtained  from 


a  smooth  curve  drawn  thru  the  points  at  a  distance  coinciding  with  the  toe  of 
the  test  slope  was  interpolated  as  the  wave  height  value  for  that  particular 
generator  setting  and  structure  slope.  Using  the  wave  height  at  the  structure 
toe,  water  depth,  and  wave  period,  deepwater  wave  height  was  computed  from 
linear  wave  theory  via  an  inverse  transformation.  Original  runup  results  were 
plotted  using  the  deepwater  wave  height  rather  than  the  measured  wave  height. 
Wave  periods  for  the  test  data  ranged  from  0.72  to  5.00  seconds  while  wave 
height  ranged  from  0.01  to  0.19  meters. 

A  listing  of  the  data  test  conditions  is  provided  in  Table  1. 

Analysis  Procedure 

As  in  all  approaches  to  evaluating  laboratory  data,  there  are  two  meth¬ 
ods  of  determining  the  important  dimensionless  groupings  of  variables  for  data 
presentation.  A  fundamental  method  for  obtaining  important  dimensionless 
groupings  of  parameters  often  utilized  where  the  physics  of  the  processes  are 
unknown  or  not  well  understood  is  the  Buckingham  Pi  method  (see  for  example 
White  (1979)).  This  method  or  a  variation  of  it  has  been  utilized  in  various 
studies  of  runup  (see  for  example  Technical  Advisory  Committee  on  Protection 
Against  Inundation  (1974)).  The  resulting  dimensionless  groupings  of  runup 
variables  for  the  case  of  normal  wave  incidence  and  linear  smooth  slope  struc¬ 
tures  of  unknown  functional  form  is  as  follows: 

R/H  =■  f[H/L0,  H/d.  pH'/oT*,  pHz/\iT,A)  U) 

where  R  -  runup,  H  -  wave  height  at  toe  of  structure,  Lg  -  deepwater  wave 
length,  d  -  water  depth  at  toe  of  structure,  p  -  density  of  fluid,  a  -  surface 
tension,  p.  -  dynamic  viscosity,  and  A  -  structure  slope.  The  third  dimension¬ 
less  grouping  on  the  RHS  (right  hand  side)  of  equation  1  is  the  Weber  Number 
for  oscillatory  flow  which  is  of  negligible  importance  for  models  of  reason¬ 
able  size.  The  fourth  dimensionless  grouping  on  the  RHS  of  Equation  1  is  the 
Reynolds  Number  for  oscillatory  flow.  Projected  effects  of  Reynolds  Number 
are  unknown  in  the  present  data  set.  Actual  scale  effect  in  runup  studies  has 
been  investigated  by  various  authors  (Saville  (1958),  Fuhrboter  (1986)),  but 
results  to  date  have  been  insufficient  to  define  such  effects  well. 
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Of  the  remaining  three  parameters  in  Equation  1,  Iribarren  (1949)  first 
noted  the  importance  of  the  combination  of  wave  steepness  (H/L*,)  and  structure 
slope  (A)  in  one  parameter,  the  Iribarren  number  (tan  (A)/(H/L0)1/2) .  Various 
researchers  (Galvin  (1972),  Battjes  (1974),  Hunt  (1959))  have  noted  the  impor¬ 
tance  of  this  parameter  in  both  the  breaking  process  and  resulting  runup  on 
beaches.  For  mild  slope  structures  Hunt  (1959)  recognized  that  the  relative 
runup  was  proportional  to  the  Iribarren  number  under  breaking  wave  conditions. 
Battjes  (1974)  provides  a  physical  explanation  for  the  relationship  between 
runup  and  the  Iribarren  number  for  the  case  of  mild  slope  structures.  A  limi¬ 
tation  to  this  equation  can  be  seen  for  the  case  of  steep  sloped  structures 
where  the  Irribarren  number  approaches  infinity.  In  the  present  paper  it  is 
desired  to  provide  a  unified  approach  for  wave  runup  on  both  steep  and  mild 
slope  structures,  therefore  a  sin(A)  term  was  used  in  place  of  the  tan(A)  of 
the  Iribarren  number  and  hence  this  modified  dimensionless  grouping  will  be 
referred  to  as  the  modified  Iribarren  number.  It  should  be  noted  that  a 
slight  refinement  of  Battjes  (1974)  arguments  (i.e.  considering  wavelength 
defined  along  the  slope  rather  than  in  the  horizontal  plane)  will  lead  to  the 
sin(A)  term  used  here.  This  modification  of  the  Iribarren  number  is  consis¬ 
tent  with  various  criteria  for  delineation  of  the  zone  between  breaking  and 
nonbreaking.  For  example,  Munk  and  Wimbush  (1969),  provide  an  expression 
based  on  linear  wave  theory  for  breaking  on  a  slope  in  which  the  downslope 
component  of  the  particle  acceleration  cannot  exceed  g  sin(A).  As  Battjes 
(1974)  noted,  with  proper  accounting  of  the  reflected  wave  height,  the  Munk 
and  Wimbush  (1969)  criteria  can  be  written  as: 


[sin  (A)  /  (H/L0) 1/2]  c  »  (2x) 1/2 


(2) 


< 

Miche  (1951)  using  linear  wave  theory  also  derived  a  kinematic  criterion  for  , 

the  limiting  conditions  of  non-breaking  on  a  plane  slope  extending  to  deep  , 

water.  His  criterion  is  given  by: 

{Hj La)  e  *  (sin2  (A)  /n)  (2A/ it) 1/2  for  Ain/4  (3) 


This  expression  can  be  reformulated  as: 
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[sin  (A)  /  (H/La) 1/2] c  -  (n/ Ka) 1/2  (x/2A) 1/4 


(4) 


where  K,  is  the  linear  shoaling  coefficient  if  one  assumes  the  expression  good 
for  all  depths.  Keller  (1961)  finds  a  similar  expression  for  limiting  condi¬ 
tions  for  non-breaking  where: 


(H/L0)c  -  JT.MV 2 ?t)  (2A/rc) 1/2 


(5) 


which  can  be  rewritten  as: 


[A/  (H/L0) l'2]  e  -  <2n/tf,)1/2  (it/2 A)1/4 


(6) 


Keller's  expression  is  based  on  a  nonlinear  shallow  water  theory  and  therefore 
might  be  considered  more  valid  than  the  expression  of  Miche.  For  small  slopes 
A  -  sin(A)  and  hence  Keller's  criteria,  Equation  6  ,  is  within  a  constant  of 
Miche's  expression,  Equation  4.  For  large  slopes  (i.e.  limiting  case 
A  -  w/2) ,  the  factor  sin(A)  differs  from  A  by  50  percent. 

An  additional  dimensionless  grouping  d/Lo  can  be  formed  by  dividing  the 
first  RHS  parameter  grouping  by  the  second  parameter  grouping.  This  grouping 
of  parameters  has  the  advantage  of  delineating  the  relative  water  depth  which 
the  structure  is  situated  in.  An  important  parameter  for  free  surface  flows 
not  explicitly  mentioned  above  is  Froude  number,  which  for  oscillatory  flows 
in  deep  water  can  be  represented  as  the  multiplication  of  the  first  two  RHS 
groupings  in  Equation  1. 

A  second  method  of  obtaining  dimensionless  variable  groupings  of  impor¬ 
tance  is  by  casting  the  physical  equations  into  dimensionless  form.  The  basic 
equations  of  fluid  dynamics  would  point  out  the  importance  of  Reynolds  number, 
Froude  number,  and  Webber  number  as  before.  To  obtain  further  groupings  of 
importance,  a  direct  look  at  physical  equations  for  runup  due  to  nonbreaking 
waves  is  called  for. 

Koh  and  LeMehaute  (1966)  have  suggested  a  runup  equation  for  sloped 
structures  of  the  form: 
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R/W  -  (n/2 A)1/2  +  {kH/ L)  (1/ tanh(lcd) )  (l+3/4sinh2  (kd)  -1/4 cosh2  (kd) ) 


(7) 


The  first  term  is  based  on  an  earlier  linear  expression  derived  by  Miche 
(1951)  for  deep  water  conditions,  while  the  following  terms  are  based  on 
Miche 's  (1951)  approximation  to  non-breaking  runup  on  vertical  walls  for  non¬ 
linear  wave  theory.  Except  for  a  missing  linear  shoaling  coefficient  the 
first  part  of  this  expression  agrees  with  chat  of  Keller  (1961)  which  was 
derived  for  non-plane  beds  with  non  steep  slopes. 

Keller  and  Keller  (1965)  derived  an  expression  for  the  case  of  a  plane 
slope  and  horizontal  bottom  using  linear  long  wave  theory  with  the  result: 

R/H ‘[JO  ( (kad) 1/2 /A)  *  Ji  ((k0d)1/2/A)]-1/2  (8) 


where  JO  and  J1  are  Bessel  functions  of  the  zero  and  first  order  ’-espectively . 

In  both  the  above  expressions  for  non-breaking  wave  runup,  the  relative 
runup  is  seen  to  be  of  the  functional  form: 

R/H  =»  f  (  A,  d/La,  H/L0)  (9) 

The  grouping  of  the  parameters  on  the  RHS  of  equation  9  can  also  be  combined 
and  reexpressed  as  before: 


R/H  =  f  (A/ (H/La) 1/2 ,  d/La) 


(10) 


where  the  first  dimensionless  grouping  is  of  a  similar  form  to  the  modified 
Iribarren  number  used  previously. 

In  an  attempt  to  unify  the  presentation  of  relative  runup  on  smooth 
slopes,  the  primary  independent  variable  of  importance  was  chosen  to  be  the 
modified  Iribarren  number. 


Results 


Relative  runup  plots  for  nine  slopes  ranging  from  vertical  to  1  on  10 
have  been  presented  in  Figures  1  thru  9.  The  expected  criteria  for  breaking 
waves  given  by  the  modified  Hunt  (1959)  expression 


R/H  -  sin  (A)  /  (H/ Lc) 1/2 


(ID 


is  superimposed  on  the  plots  along  with  an  upper  limit  expression  of  the 
Miche-Keller  form  (Equations  2  and  4)  for  non-breaking  wave  limit.  The 
expression  for  this  upper  limit  found  to  be  most  consistent  with  the  data  is 
given  by: 


R/H  =  (2ti)  1/2  (it/2A) 1/4 

If  the  shoaling  coefficient  Ka  is  assumed  to  be  unity,  this  expression  is 
consistent  with  Eq.  4  (except  for  a  constant  (-  2.0))  and  with  Eq.  6  (using 
Sin(A)  rather  than  A).  Table  2  presents  this  relative  runup  upper  limit  ver¬ 
sus  structure  slope  for  the  slopes  investigated  in  this  study.  The  rational¬ 
ization  for  this  approach  to  maximum  wave  runup  is  that  within  the  realm  of 
breaking  waves  (on  the  slope)  the  relative  runup  should  follow  the  modified 
Iribarren  number.  As  the  modified  Iribarren  number  is  increased  and  enters 
the  zone  of  non-breaking,  the  relative  runup  should  decrease,  therefore  limits 
provided  should  envelope  runup  data. 

In  all  cases  the  data  in  the  breaking  zone  portion  of  the  curve  follow 
the  modified  Iribarren  number  relationship  well.  Additionally,  the  critical 
transition  point  for  nonbreaking  describes  well  the  upper  limit  of  the  rela¬ 
tive  runup  except  for  the  1  on  10  and  vertical  slope  where  it  overestimates. 
The  reason  for  this  overestimation  on  the  1  on  10  slope  is  not  entirely  known 
but  believed  to  be  due  in  part  to  the  greater  effect  of  friction  as  slope  gets 
milder  and  the  consequent  opportunity  for  a  viscous  boundary  layer  to  develop. 

The  upper  limit  for  relative  runup  on  a  smooth  vertical  slope  according 
to  Eq.  (12)  is: 
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R/H  =  (2r)1/j  =  2.5 


<13) 


A  simple  periodic  standing  wave  on  a  vertical  slope  would  produce  an  expres¬ 
sion  with  relative  runup  R/H  -  1.0.  Additionally,  Wallace  (1964)  provided  a 
method  to  numerically  calculate  the  runup  for  solitary  waves  (which  can  be 
considered  to  be  a  limiting  case  ef  periodic  waves  in  shallow  water) .  For 
vertical  walls  with  H/d  >0.15  his  method  produces  an  estimated  R/H  -  2.5  in 
accord  with  the  above  proposed  criterion.  A  value  of  relative  runup  R/H  =2.5 
on  vertical  waves  is  not  inconsistent  with  laboratory  findings  of  Takada 
(1974)  or  the  original  recommended  value  of  R/H  **  3.0  proposed  by  Hunt  (1959) 
for  surging  (non-breaking)  waves. 

Although  graphs  have  not  been  included  in  the  present  paper,  the 
approach  has  been  used  on  one  set  of  data  with  a  1  on  3  slope  in  29.5  cm  water 
depth  and  is  consistent  with  the  results  provided  herein. 

i 

Conclusions 

A  method  is  presented  for  estimating  the  upper  limit  of  periodic  wave 
runup  on  smooth  sloped  coastal  structures.  The  method  is  simple  and  can  be 

I 

applied  in  an  unambiguous  manner  to  a  wide  range  of  slopes.  Essentially,  the 
method  is  composed  of  two  elements,  a  modified  surf  parameter  and  a  limiting 

i 

value  for  relative  runup  caused  by  non-breaking  wave  conditions.  The  modified 
surf  parameter  allows  vertical  structures  to  be  included  in  the  analysis  which 
would  be  impossible  with  the  standard  surf  parameter,  Battjes  (1974).  For 
structures  with  steep  slopes  or  vertical,  where  non-breaking  wave  conditions 
dominate,  the  method  provides  a  logical  envelope  for  the  runup  data.  On  some¬ 
what  flatter  slopes  where  a  mixture  of  wave  conditions  occur,  the  method 

f 

predicts  the  runup  of  breaking  waves  quite  well  and  gives  a  reasonably  conser¬ 
vative  envelope  for  the  runup  of  non-breaking  waves.  Minor  modifications  of 
the  method  might  be  required  for  non-breaking  conditions  on  flat  slopes  to 
ensure  against  gross  over-prediction. 


< 
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Table  1 

Summary  of  Test  Conditions 


Structure 

Slone 

Wave 
Period 
( sec') 

R 

(cm’) 

H/d 

d/L 

0 

Number  of 
Observations 

vertical 

0.9-12.2 

0.72-4.71 

1.3-15.5 

.024-. 321 

0.11- .471 

26 

1  on  0.5 

0.9-18.6 

0.72-4.71 

1.1-37.2 

.025- .488 

.011- .471 

33 

1  on  1.0 

0.9-18.6 

0.72-4.71 

0.9-42.8 

.022-. 488 

k 

.011-. 471 

32 

1  on  1.5 

1.2-17.1 

0.72-4.71 

1.2-39.6 

.030- .448 

.011- .658 

45 

1  on  2.25 

0.9-17.4 

0.72-4.71 

1.5-45.8 

.028-. 456 

.011- .658 

51 

1  on  3 

0.9-17.7 

0.72-4.71 

2.4-48.4 

.023-. 464 

.011- .658 

49 

1  on  4 

0.9-18.3 

0.72-4.71 

2.3-47.5 

.028- .480 

.011- .658 

51 

1  on  6 

0.8-17.8 

0.72-4.71 

1.8-37.2 

.021- .468 

.011- .658 

51 

1  on  10 

0.5-12.1 

0.72-4.71 

0.8-8. 9 

.013-. 480 

.011- .471 

28 

Table  2 

Maximum  Relative  Runup  Versus  Structure  Slope 


Slope 


R/H  (equation  12) 


> 


vertical 

2.50 

1:0.5 

2.73 

1:1 

2.97 

1:1.5 

3.19 

1:2.25 

3.49 

1:3 

3.70 

1:4 

3.95 

1:6 

4.39 

1:10 

5.00 

Slope  1:10  (d/Lo  =  0. 01  1 -0.471 ) 


Slope  1 :4 


(d/Lo  =  0.0  1  1  -0.658) 


6  20 

sinA/ sqrt(H/Lo) 

Figure  3,  Relative  runup  versus  modified 
Iribarren  number,  slope  1:4 

Slope  1:3  (d/Lo  =  0. 01  1 -0.658) 
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Slope  1:2.25 


(d/Lo  =  0.01  1  -0.658) 


Figure  5.  Relative  runup  versus  modified 
Iribarren  number,  slope  1:2.25 


sinA/ sqrt(H/Lo) 

Figure  6.  Relative  runup  versus  modified 
Iribarren  number,  slope  1:1.5 


Slope  1 : 1 


(d/Lo  =  O.Q  1  1  -0.471  ) 


H/d=.022-.488 


Figure  7.  Relative  runup  versus  modified 
Iribarren  number,  slope  1:1 

Slope  1:0.5  (d/Lo  =  0.0 1  1  -0.47 1 ) 


H/d  =.025-. 488 


sinA/ sqrt(H/Lo) 

Figure  8.  Relative  runup  versus  modified 
Iribarren  number,  slope  1:0.5 


R/H 


Slope  Vertical  (d/Lo  =  0. 01  1 —0.4-71 ) 


Figure  9.  Relative  runup  versus  modified 
Iribarren  number,  slope  vertical 


Appendix  A:  Notation 


The  following  symbols  are  used  in  this  paper: 

-  angle  of  structure  slope  with  horizontal; 

-  wave  height  (at  toe  of  structure),  deepwater  wave  height; 

-  linear  wave  theory  shoaling  coefficient; 

-  wavelength,  deepwater  wavelength; 

-  runup; 

-  wave  period; 

-  water  depth  (at  toe  of  structure) 

-  function  of; 

-  gravitational  acceleration; 

-  wavenumber;  deep  water  wavenumber 

-  water  density; 

-  dynamic  viscosity;  and 

-  surface  tension. 


